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Abstract 

We study SU{N) S{U{N/2) x U{N/2)) symmetry breaking in AA = 4 
SYM via AdS/CFT correspondence. Two stacks of ^ parallel D3 branes 
are separated by a distance 2d in the transverse directions. This leads to 
an interaction between different l-wave dilatonic KK modes when expanded 
in harmonics. We calculate certain two point correlation functions of the 
dual gauge theory in the decoupling limit ^ 1 . Due to mode mixing, 
the diagonal correlation functions have jj conformal-like correction as well 
as deformation terms. The off-diagonal correlators are also nonvanishing and 
their leading order is jj. We discuss briefly the spectrum of the glueball 
exitations. 



1 Introduction 

A recent development in the understanding of the string theory was a conjecture by 
Maldacena according to which, in some appropriate low energy limit, there is a 
duality between A/" = 4 supersymmetric SU{N) gauge theory and type II strings on 
AdSd+i X JH background. This is nothing but a manifestation of the old ideas 0,0 
involving the deep connection between gauge theories and strings. Therefore, from 
this point of view it is clear that many problems in large N gauge theories can be 
studied via supergravity in AdS spaces|. The precise meaning and the explicit rules 

-'^e-mails: rrachkov@sfu.ca; rash@phys.um-sofia.bg, on leave of absence from Dept. of Physics, 
Sofia University, 1164 Sofia, Bulgaria 
^e-mail: kviswana@sfu.ca 
•^For a recent review see for instance [p|. 
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for relating AdS supergravity and the type IIB string theory were given in and . 
Namely, the explicit form of the correspondence between the correlation functions 
of the operators in conformal field theory (CFT) living on the asymptotic conformal 
boundary of the d + 1 dimensional anti-de Sitter space and type II supergravity 
action evaluated on the classical solutions was proposed to be: 



{exp j %0,)cFT = exp {-S[^']) 



where OiS are conformal operators in CFT and $9 = ^i^^^^^ is the boundary value 
of the classical solutions which serve as a source for the primary operators in CFT. 
The proposed correspondence has been examined in a number of calculations of the 
correlation functions of operators in CFT using the classical solutions of various 
fields from the spectrum of AdS supergravity. These, even incomplete proofs of the 
validity of the holographic principle, supply us with a powerful tool for studying 
strongly coupled gauge theories. Most of the intensive investigations of AdS/CFT 
correspondence however are for the case of supergravity in pure AdS background, 
i.e. at the origin of the moduli space where there is an enhanced superconformal 
symmetry and many problems can be managed exactly even at strong coupling. 

Departure from the origin of the moduli space means allowing scalar fields to have 
non-zero expectation values and therefore to move to the Coulomb branch of the 
theory. The vacua in the Coulomb branch are maximally supersymmetric but some 
of the gauge symmetry is broken by the expectation value of scalar fields. These cases 
correspond to certain solutions of the supergravity field equations in which the bulk 
space geometry approaches AdS near the boundary, but differs from it in the interior. 
From holographic correspondence point of view such geometries are related either to 
a relevant deformation of the conformal field theories, or to existence of states of the 
CFT with vacuum expectation value which sets the energy scale. In a recent paper 
by Klebanov and Witten a more general investigation of the symmetry breaking 
is given and apphed to several examples. There are recent papers studying the 



Coulomb branch of the gauge theories in four dimensions |jT^, |T0], [TT], ^ |T6[ which 
are based on continuous spherically symmetric distributions of large number of Wi 
branes. They, in fact, extend the results of [12, 1^ by investigating different probes 



in the bulk of the AdS space and their " holographic images" on the boundary. Since 
the effective background is asymptotically AdS;^ x one can expand the fields in 
spherical harmonics and solve the supergravity equations of motion for the Kaluza- 
Klein (KK) modes. In this case, due to spherical symmetry of the distribution of 
the D?> branes the different KK modes decouple, thus enabling one to find exact 
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solutions of the field equations using AdS/CFT correspondence to obtain the two 
point correlation functions. 

In the present paper we will consider the simplest example of two stacks of equal 
number coincident D3 branes separated by a distance 2d, which corresponds to a 
vacuum state of the gauge theory where the gauge SU {N) symmetry is broken down 
to S{U{N/2) X U{N/2)) by the expectation value of scalar fields. 

In Section 2 we consider the supergravity solutions in the presence of two center 
D3 brane geometry and the reduction of the 10 dimensional dilaton action to five 
dimensional one by integrating over S^. It turns out that the equations of motion 
for the dilatonic KK modes are coupled. This coupling modifies the two-point 
correlation functions, as discussed, in Section 3. The diagonal correlation functions 
of the operators in the gauge theory are modified by a correction term to the 
conformal-like part as a consequence of the coupling between the different modes. 
In contrast to the conformal case and the case of uniformly distributed D3-branes, 
the correlation functions {0"'{x)0^~^'^{y)) are nonvanishing and are of order 1/N 
in our approximation. In addition, all correlators have a deformation term which, 
after a Wick rotation, exhibits a pole structure. We interpret the appearance of 
double poles as a degenerate spectrum of glueballs. All these questions are subject 
of consideration in Section 4. We summarize, interpret and comment on our results 
in the Concluding section. 



2 Two stacks of D3 branes: equations and solu- 
tions 

Let us consider the case of two center D3 branes . In this large number 

of D3 branes are separated into two parallel stacks of N/2 branes each. Their 
separation distance introduces an energy scale in the theory. The SU{N) gauge 
symmetry is broken down to S{U{N/2) x U{N/2)). iFiom the gauge theory side 
this corresponds to Higgsing of the gauge theory by scalar fields: 



X = "^-^ y (1) 




where / is the {N/2) x {N/2) identity matrix and di is the position of the i-th stack 
in the transverse direction. Let us parametrize the moduli space in terms of gauge 
invariant fields as in |jl5|, ||: 

= Cl,,Tr {X^^ . . . X^^) (2) 
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where C/^ is a totally symmetric traceless rank k, SO (6) tensor and X*'s are as 
in (|l]). The Higgsing gives expectation values for the chiral operators which 



can be expressed in terms of order parameters di H, 15 



{O') ~ n + n (^2), - traces 



k=l 1=1 



To describe correlation functions on the Coulomb branch let us start with the 
ten dimensional Type II supergravity in the two-centered D3 brane background: 
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ds^ = H 5 [-dt + dxi + dx2 + ^0:3) + H^^ dy^ (3) 
where H is the Green's function with two separated sources of equal charge: 

(4) 




In y are the coordinates normal to the branes 0. 

Let us discuss briefly the geometry of this configuration and the decoupling limit. 
Each of the two stacks develops a throat (geometry of which is AdS) but they are 
located at different points of 5*^, separated by a distance 2d. In each throat the 
dilaton field can be expanded in 5*^ scalar spherical harmonics which will have non 
vanishing overlap with the corresponding expansion around the other throat. For 
small d (compared to the radius of AdS) and large r the two throats meet and in the 
limit r ^ 00 share a flat four dimensional boundary. In this case it is natural to make 
an appropriate shift in the, ^-direction, choosing the origin so that di = —d2 = d 
and therefore: 

r4 / 1 1 \ 

(5) 

which allows us to consider the problem as effectively single centered. Due to inter- 
action between the different KK modes in the bulk, it is expected that their images 
on the boundary result in a mixing between the correlation functions of operators 




^We call these operators chiral in the sense of [|l5| 

^It is clear that at large r the space is asymptotically AdS^ x S^. Since we consider the large 



N case we have neglected the one in H. L = Airgga' N. 
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of different conformal dimensions. We would like to point out that starting from 
Type TIB string theory and taking the decoupling limit a' ^ while keeping d/a' 
fixed, we obtain an A/" = 4, 4-dimensional SYM theory with gauge group SU (N) 
spontaneously broken down to S{U{N/2) x U{N/2)) by Higgs scalar expectation 
value in the adjoint representation of SU{N). The Higgs scalars parametrize the 
positions of the stacks. We use the supergravity approximation which is valid for 
d « L^. In general one must distinguish between three different regions: the 
first one is near the location of the branes. Here the throats are well separated. 
The second is where the throats come together, but this region is still far from the 
boundary. The last region is near the boundary where the two throats effectively 
coincide. Since our coordinate system is symmetric with respect to the positions 
of the stacks, the expansion over spherical harmonics qualitatively changes on the 
boundary r = \d\. There we will apply the matching condition technique in order to 
ensure well defined smooth solutions in the bulk. Therefore, using this choice of the 
origin of the coordinate system one can expand the harmonic function H in Taylor 
series in the form: 



where = ^ Cil^.i^iV^^ ■ ■ ■ V^^'' larger or smaller of r = \y\ and d 

respectively. 

Let us now consider the kinetic part of the lOD dilaton action: 







After integration by parts the above action takes the form: 



t$ _ 

10 — 



Qkm I Q 
-^10 + 



(8) 



where: 





AdS 



and 




(10) 



dAdS 
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One can substitute for the dilaton field its expansion in 5*^ scalar spherical harmonics 



=^<l>^(f,r)y^(fi5) (11) 

7=0 

where r = \y\, and to obtain the five dimensional action in the background metric (^) 
we integrate over 5*^ . As it was noted in a number of papers (see for instance ||T7|), 
we have an infinite tower of KK modes and there are no fundamental physical 
reasons to truncate. The resulting boundary action takes the form: 




(14) 



A 



ij 



L 



1=1 ^ ^ 



r) 2^-i(/+l)(/ + 2) 

21 

a{2l,I,J){C^'C'C'') 



As in |]T3[ we use the following normalization: 

1 kL\ki\kj\ 



aim, J) 



(S + 2)!2^-i'ai!a2!«3! 



S = -(/cl + /c/ + /cj); «! = S - /cl; a2 = T, - kj; as = S 



ki 



(15) 



(16) 



We will see that the last term in ([121) produces only local terms on the boundary 
and therefore is irrelevant for our considerations. 



^We will use the notations and normalization of the spherical harmonics as in |I5| 



6 



The kinetic term can be easily reduced to five dimensional effective action in the 
form [H (from now on we will denote by $): 



S. 



kin 



d xdr < $ 



dQ^H \ Vl + Vr,—drr^dr 



+ ( I dn^Y^Y-^ ^^^drr^dr^^ + 4 I /" dn^Y^V%Y-^] ^^^A (17) 



It is straightforward to derive the equations of motion following from (|T^ in the 
form: 



'1^ 



dn^Y^Y^ A,$-^+— / dn^Y^wl.Y-^ + 



J dn^HY^Y-^^ Vp^ + (^j dn^HY^^ Vj$ = 0, (19) 

where = {l/r^)dr{r^dr). 

Since the metric has different forms in the regions ofr < d and r > dwe will have 
two systems of equations - one for each region. Taking into account the expression 
for H (^) and the orthogonality of the spherical harmonics (with normalization as 
in [|T3|) the equations of motion are given by: 

a) r > d region 



1=1 ^ ^ 



k). 



(20) 



A,$"(r, k) - ^^^^y^<l>"(r, k) - ^A;2$"(r, k) 



1/ 

»,4 



A;22"-i(n + l)(n + 2)^22'(2/ + l) (- j ^ ^2/,n,m^ 

1=1 ^'"^ m=\2l-n\ 



{r,k) 



21 



+ ^k^2^i^2l + l) (-) 5''^'''^r,k). (21) 
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h) r < d region 



A,,<|.(r, fc) - -e^r, k) = ^k'Y,j— (-) $^^(r, fc) (22) 



21 



^A;'2"-i(n + l)(n + 2)^22'(2/ + l) (^^) ^ ^^''"•'"^'"(r, A?) 

i=l m=\2l~n\ 

+ _fc222'(2/ + l) (^-j 52Z,n^(r,fc) (23) 



In the above expressions we used the following notation: 

A^'^'-^ = a(L, J, J) (C^C^C^) (24) 

and a(L, /, J) and < C^C^C"^ > are as in eq.(|T6|). We see that the equation for the 
different KK modes of the dilaton field are coupled. In general, according to the 
AdS/CFT correspondence, the operators are mapped to conformal fluctuations in 
the metrics of AdS^ and S^. In our case this implies that the dilaton action projected 
on the boundary will produce correlation functions of the operators ~ tr [F'^X^). 
In contrast to the uniformly distributed D3 branes case, there is coupling between 
the different KK modes, which signals a nontrivial contribution to the corresponding 
correlations functions. The mixing between the KK modes means that, due to the 
fact that we have a broken conformal symmetry, the sources coming from gravity 
solutions will produce in general mixed correlation functions. We will see how it 
happens in the next section. 

The field equations in the two regions can be transformed into Bessel equations 
by the following substitutions: 

a) for r > d let us define: 

r 

i = = ^k'; (e < 1) 
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then the equations ^) become: 

^ ^ / — 1 



i=l m=\2l—n\ 



+ 



22'(2/ + 1) 



21 



b) For r < d, defining: 



2 _ 7, 2 



(^<1) 



the equations become: 



1 / /I \ CXD ^ 

/" (^) + -J'iv) - (i + 7^) /(^) = E (7TI)^^"^-(^) 

^ / — 1 



(26) 



(27) 



/:(^) + -m - f 1 + /n(^) 

= 2"-^(n + l)(n + 2)E ^j^r E ^''•"''"/-(^) 

i=l m=|2i— n| 

^ 2^H2/ + 1) .,,^,,„ (28) 

where primes on / denote differentiation with respect to its argument. The terms on 
the right hand side of the equations above involve couphng to KK modes. For small d 
and large A^, we can solve these equations perturbatively. One can see that the right 
hand side of both systems of equations is readily expanded in d^ / = n'^m^jj, 
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which in the decouphng hmit is small. Here, niw = is the W boson mass, 
generated by symmetry breaking. This is precisely the region of validity of our 
approximation, i.e. for energies below the Pi^-boson mass. In the zeroth order we 
ignore the coupling of the modes and these solutions are used as inhomogeneous 
terms in obtaining the first order solution. 

Let us consider the zeroth order approximation. The solutions of the homoge- 
neous system for r > d are: 



fniik) 



n = 0, 1,... 



4+2(0 

and therefore the general solution for the dilatonic modes is: 

= e [i^„+2«)/if (^) + /„+2«)/xi°n^) 

where fin\k) and iJn\k) are yet to be determined coefficients. 
The solutions in r < d region are: 



(29a) 



(30) 



fniV^k) 



n = 0,1, 



(31) 



or, taking into account that rj = 1/^: 



ir.+2(|)«i°nfc)+4+2(|)«r(^) 



(32) 



In general, our solution must be regular in the interior of AdS space so that the 
solutions for r > d must be extended in r < d region to a regular solution. The only 
regular solution in this region is: 



(33) 



In addition, the two solutions ( PPD and (^) must be matched on the boundary of 
the two regions (^ = 1): 



(34) 
(35) 
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These conditions uniquely determine the coefficients jjih^ and 



24+2(«:)/;+2('^) 



^«+2(/t)/;+2(^) - -^n+2(^)-^n+2(^) 



(36) 



(37) 



To determine the last unknown coefficient jln \ note that the scalar field near the 
boundary behaves as 



where Qb{k) is the boundary value of Q and u is the spectral parameter in the 
equations of motion. Applying the above formula to the n-th dilatonic KK mode 
and using the behaviour of the Bessel functions of small argument, one can see that 
the only contribution to the above behaviour comes from Kn+2 and therefore: 



or: 



{n + l)\ ^K^~(n+2) ^^^^^ 



n+2 



n{k) 



(38) 



(n + 1)! V2- 

With this, all the coefficients in the zeroth order solutions are completely determined. 
We subsitute the zeroth order solutions on the right hand side of ( p5| , p6D and ( p7| , p8D 
to find the first order correction. The equations for the first order (/ = 1) are as 
follows: 

a) r > d region 



1+ 



f{i) 



1 



1 

/n(0 



m)f^f\k)+h{Of^T{k) 



(0), 



(39) 



^2,n,n-2 (^k£) f^^^l^^k) + /„ (O/i^l^ (^) 
+ A''-'- (K^+2iM\k) + In+2{i)l^T{k) 

l2ik)+In+d)l^% 



^^2,n,n+2 ( R^^,^^) fl^^l^{k) + In+d)l&2{k) 



(40) 
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h) r < d region 



fiv) + If'iv) 

V 



2k^ 



m) + \m) - f 1 + ) uv) 



^ 2"+i(n+l)(n + 2) .2 



(0) 



Solving these equations we find tliat the general solutions for r > d are[] 



2,n,n 



(0) . t3i^ 

3r /tA , Q^^ , t\c2t \ ,,(0) 



e%,_i(o + 3(n+ir/„(e)j/x, 

- + 3(n + l)f ir„(0) 



For the region r < d, the regular solutions are: 



, 2"+^(n + l)(n + 2) 



A 



2,n,n-2 -3 



(0) 



+^2'"'" (f /„_i(r)) + 3(n + l)f /„(r/)) + A^'^'^+^f /„+5(^)4' 



(0) 
+2 



^See Appendix for details. 
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Once again we will use the required behaviour of the solutions near the boundary. 
Expanding the Bessel functions in (^3]) and (|4^) for small argument we find that: 



[n 



1)! 



{k) + 2 



n+l 



[n + r 



{k) 



(47) 



/}«(«:) = 2 



As for the zeroth order solutions, the coefficients fin\K) are determined by the 
boundary data. There is, however, an important difference - these coefficients are 
determined by a linear combination of ^h{k) and $^'''^(fc). The solutions for n-th 
dilatonic KK mode contain two as yet undetermined coefficients, but they can be 
determined by the matching conditions on the boundary of the two regions. We 
will see later that these coefficients will play an important role in interpreting the 
correlation functions and therefore on the physics of the Coulomb branch. 

To conclude this section let us summarize the results we have found. 

Starting from lOD dilaton action in 2 stacks D3 brane background we derived 
the coupled equations of motion for the KK harmonics. We have obtained explicit 
solutions to first order in ■\. The arbitrary coefficients in the solutions were deter- 
mined by applying the boundary conditions at r = d as well as from the appropriate 
asymptotic properties. The most interesting results are contained in eqns(47) and 
(48) which reflect the coupling of the n and n + 2, KK modes. This leads to non 
conformal correction terms to the correlation functions. 



3 Correlation functions 

As mentioned in the Introduction two stacks of parallel, separated D3 branes pre- 
serve supersymmetry, since the Poincare supersymmetries of the gauge theory are 
maintained but super conformal invariance is broken by the Higgsing. In the case 
under consideration, the SU{N) symmetry is broken down to S{U{N/2) x U{N/2)). 
This will modify the correlation functions which we study in this section. 

lYiom. the gauge theory side, according to AdS/CFT correspondence, evaluation 
of the action on the classical solutions will produce the correlation functions of the 
operators tri^F'^X^^ (we use the notations of 0: = C/^ j^X*^ . . . X*''). As 



shown in the Introduction, the nontrivial contributions will come from the terms of 
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the following formQ : 



Sb = -- lim [ d^xe-^^^d,^l^ 

2 e-*0 J ^ 



(49) 



We saw in the previous Section that there are two undetermined coefficients 
fin\K,) and an\K,) for every n. In order to proceed with the derivation of the 
correlation functions in the gauge theory we ffist determine these. For this purpose, 
we have to impose on the first order solutions the matching conditions (|34| , ^5D on 
the boundary of the two regions = 1) (note that these are imposed now on the 
ffist order solutions). After tedious but straightforward calculations we find the 
following expression for yUn "*(/«): 



where: 



/ii'H'^) = f3n,n^2{K)^r'0) + f3nA^^)^b fi) + /9„,„+2 («:) (fc) 



(50) 



(51) 



+ 2"(n + l)(n + 2)^2'"'" [{kK^^i + 3(n + l)Kn) l'„ 



+2 



X 



(n+l)!4+2/;+2 



(52) 



2"(n + l)(n + 2)/l2'"''^ 
" (n + l)\ 



[kL^^i + 3(n + 1)4) /; 



+2 



+ (34_i + + ^^^^4 + 3(n + 1)/:^ 



t2 jt •! 
^n+2^n+2 



(53) 



/?n,n+2(«:)=C+2('^)+/5S+2('^) 



(54) 



^One can easily check that the terms which survive in the hmit ^ — » coming from Ai^'^d^^^ 
contain only positive integer powers of k and hence are contact terms. 
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(1) 

n+2 



2"(n + l)(n + 2)A2'"'^+2 



(n + 1)! 

In+2 

(n + 2)(n + 3) 
2'*(n + l)(n + 2)A2'"'"+2 



2(n + 4) (i^n+2/;+2 + ^;+24+2) (2 



1 



In+2ln+2 



(55) 



/-'n,n+2 



2(n + 3)! 



\Kln+bl'n+2^n+Aln+i 



+ (3/„+5 + /^-^n+s) In+2K'^+Jn+A\ J 7) 7 

^n+2^n+2^n+4^n+4 



(2) 



(56) 
(57) 



(1) _ 2> + l)(n + 2)A^ 



2,n,n-2 



X 



X 



-'n+2-'n+2 



(58) 



(2) _ 2>+l)(n + 2)^W 

A^n,n-2 — 2(n-l)! L'^-'n-l-'n+2-^n-'n 

+ (34_2 + /n+2i^;/n] , , „ 

-'n+2-'n+2-'n-'n 



(59) 



Having completely determined these coefficients, we simply have to subsitute the 
explicit form of our first order solutions in the boundary term action and take the 
limit e — 0. To see which terms survive this limit, we must use the expansion of 
the Bessel functions for small argument. After some lengthy but straightforward 
calculations we find the following expression for the diagonal correlation functions 
in momentum space ( we have omitted terms of even integer powers of k since they 
produce in the position space (5-function and its derivatives): 

23 //t\2n+4 



+ {-iy'2 



-ir- 



(n 



+ l)!(n + 2)! \2J 



- 



„2n+i4(n + l)2(n + 2)(n + 3) ^2,n,n 

^_2^n22n+2 



In k 



\ 2n+2 



(n+l)!2 

n + 2)(5n2 + 13n + 12) . ~ ../K.\2n 



Ink 



(n- l)!(n-2)! 



(60) 
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Since we have mixing between the different KK dilatonic modes, the off-diagonal 
correlation functions are also nonvanishing. We find for them the following expres- 
sion: 



=-5{k + k') 



l^no2n+4 (n + 1) (n + 2)^ (n + 2) (n + 3) (n + 4) 



+ (-1) 



n)\ 



X 



^^2,n,n^2,n,n+2 ^ ^ L'„,„+2(«) 



(61) 



In addition to the conformal-like parts (integer powers of k times In k) the correlation 
functions contain deformations coming from the coefficient functions fin'' entering 
the solutions. Substituting the expressions for jjin^ into the boundary action and 
keeping the terms that survive ^ — > limit one can obtain the deformation functions 
Dn To simplify the final expressions we introduce the following notations: 



an(«) 



D! (2) 



(n + 1) 

,n+2(" + 2)("+4) 



^2,n,n+2 



-f^n+2-^n+2 + 



n+2-' n+2 



2-^n+2-^n+2 



In these notations, we find for D„,rn the expressions: 

2 /K\«+2 



n + 2 



R f-Y 

^nPn^n ^ 2 / 



- 2"-2(n - l)(n + l)(n + 2)(5n + 2)^2'"-2'"an/3„ (-) /5„_2,n 

2)(5n + 3)A2'"'-a^/3„ (|) + 6n-2/?n-2,n (f) 



+ 2"+i(n + l)2(n + 



-^N n+2 



(62) 



n,n+2 



n + 2 



'n,n+2 



n+4 



+ &n/?n 



n+2 



- 2"(n + l)(n + 3)(n + 4)(5n + 12)A2'"'"+2a„+2/3n+2/3n,n 
+ 2"+i(n + lf{n + 2)(5n + 3)^2,n,n a„+2fcn/3n 



22"+2(^ + lf{n + 2)=^(n + 3)(n + 4)(5n + 3)A2,n,n^2,n,n+2 I - I a„a„+2/3n (63) 
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Fourier transforming to position space we obtain the following for the two point 
correlation functions: 



(0"(f)0"(?7)) 4(n + 2)(n + 3) 



(^)2n+4 ^2 '\x-y\ 

+ 2 



2n+8 



L2 



+ br,,n{^\x-y\) + 0{l/N^) (64) 



_,„+3 c^^ (n + l)2(n + 2)^(n + 3)(n + 4) 
(^)2n+6 47r^3a''A^' tt^ 

X ^'""^'^3^ +0(VA^^) (65) 

In above, we have used = 4:7rgsa''^N, and the deformation of the conformal 
part Dn,m{jji\x — y\) is the Fourier transform of Dn,m{K)- It is well known that 
^ tr{F'^X''). In contrast to the case of spherically symmetric, continuously 
distributed £)3-branes, the mixed correlators are nonvanishing. Keeping the first 
order terms in the solutions we have found that the correlation functions contain 
corrections. The leading term in the mixed correlators, which comes from the 
first order solutions, is proportional to One can see directly from the explicit 

form of the correlation functions that in the conformal limit d ^ (or equivalently 
< X >— > 0) the correction term vanishes (the deformations Dn,m also vanish). In 
that limit the mixed correlation functions also vanish and the conformal symmetry 
will be restored. 

Let us recapitulate the results found in this Section. The physical picture that 
emerges in our study of separated two stack D3 brane background is qualitatively 
different from the case of continuously distributed D?, branes. As expected, due 
to the broken conformal symmetry, the correlation functions of the boundary SYM 
theory in its Coulomb phase have terms which vanish only in the conformal limit. 
There is a deformation term much like in the case of spherical D3 brane shell and 
conformal like correction which is of order This correction does not appear 
in the case of uniformly distributed D3 branes. Futhermore, due to interaction 
between KK modes, off diagonal correlators are nonvanishing and are of order j^. 
All correction terms vanish in the conformal limit. The term proportional to is 
not a string loop correction as we have used the classical supergravity limit. 
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4 Deformation terms and the spectrum 



In the previous section it was shown that the correlation functions on Coulomb 
branch have deformation terms. Schematically, they have the following form: 



where the coefficients q can be read off from the explicit form of Dn^m and as before: 

In the above formulae we have separated (3n,m into f3n)n, the part has in the denom- 
inator a product of a Bessel function and it derivative, and l3n,m which contains in 
the denominator the square of a Bessel function times its derivative. We note also 
that Pn contains in the denominator /„+2-^4+2- 

In order to analyze the structure of -D„,m we perform a Wick rotation k —ik. 
One can use the transformation properties of the modified Bessel functions: 

I,(-ik) ^ e-'""^ J,(k) 

K,{-ik) = ye^'^^ [J,{k)+iN,{k)] 

to obtain the relevant expression for the case of Minkowski space time. Since the 
expressions are rather complicated we will analyze only the denominators. From the 
above formulae, one can see that there are two kinds of poles. Simple poles come 
from the zeroes of: 

J„+2(k) J^+sl'^) = 0; or J ri+2{l^) J J n{l^)J'n{l^) = (66) 

and poles of order two are determined by 

[j„+2(«)j;+2(«)]' = o. 

We interpret these as describing the spectrum of glueballs with double degener- 
acy. It is reasonable to expect that the two stacks might introduce some degeneracy 
in the spectrum since they contribute identically to the boundary. The spacing 
between the different resonant frequences is given in units of It is interesting 
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to note that we obtain the same resonant frequencies as in the case of D3 brane 
shelL However the resonances are sharp. One possible way to resolve this puzzle, 
as proposed in ||^ is to consider effect of absorption by the branes which introduces 
finite width to these resonances. It would be very interesting to study the stability 
of these resonances and their eventual decay into gauge bosons. 

We show below that the deformation terms contain logarithmic corrections to the 
two point correlators due to the breaking of conformal symmetry. We will consider 
only the diagonal term in (64) (because it is leading order in -^) Then we find, 

DoM ~ (|)'/3o,o + . . . 

where ellipses stand for contributions from the mixing. In the expression for /?o,o we 
perform a Wick rotation and we find 

J2J2 ^2 ^2 

Thus we have to evaluate the expression 



(27r)4 J ^"'"^ ' (27r)2 \x - y\ 

Changing the variable k ^ k\x — y\ one finds that 

d^ke~'Us-y)p^^^^^) = ^ f dkk'j,{k)Po,o 



d 



In order to obtain the leading order contribution, one can use the series representa- 
tion for the Bessel function N2 in /3o,o 

^ I y\ 

'KN2{K.\x-y\) = J2{k\x - y\)log h . . . . 



In above the dots stand for terms of positive powers of k,\x — y\ which for a local 
theory is small (|x — y| is finite, but infinitesimally small). Let us consider first the 
term N2/ J2- The integral becomes 

d^ke-'^-^'^-^^pofi = .J'. ^ ^is / dkk^Ji{k) {log k\x-y\ + ...) 



(27r)4 J ^"''^ ( 27r)2 \x - y\ 

^ Since we want to demonstrate appearance of ^017 corrections we will not keep track of numerical 
factors. 
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To compare our results with known results we will restrict to energies below the mass 
of the VT-boson. In this case the integral over k is finite and the logarithmic correc- 
tion to the two point correlation function is reproduced. Analogous considerations 
hold in the case of N2/J2. We note one can repeat these steps for the non-diagonal 
part of 1^0,0 which will produce a logarithmic correction to the mixed term. In the 
above range of energies, from the gauge theory side, the physics is as follows ||18|,[18]. 
Since the conformal symmetry is broken by Higgs scalar field, we will have a renor- 
malization group fiow from the conformal point of SU{N) to IR fixed point with 
SU{N/2). If we consider the case of a D3 brane probe, the corresponding picture 
on the gravity side is as a string stretched between the D3 brane and the horizon of 
AdS and can be described by the Born-Infeld action which is a Yang-Mills theory 
plus higher derivative terms. Then d/a' will correspond to the mass of a VT-boson 
and the symmetry will be broken from SU{N) to S{U{N — 1) x U{1)). Now, one can 
proceed by adding more and more D3 branes to the single one until the resulting 
theory becomes S{U{N/2) x U{N/2)). All of this have correct gravity approximation 
when QgN >> 1, i.e. large N. In the large N limit only planar diagrams contribute. 
The logarithmic correction then has an interpretation as due to processes where a 
pair of virtual gauge particles is produced. The ultraviolet cut off is related to the 
gravity mass gap and in the strong coupling limit it leads to a condensate of gauge 
particles. It is shown in the above that our deformation term has poles in the Wick 
rotated k plane that correspond to bound states of glueballs. 

The main feature in our calculation is contained in the second term in (64) as well 
as in the non-vanishing values for the mixed correlators, i.e. correlators of conformal 
operators of different conformal weight. Note that these terms are proportional to 
~ J^- integrate out all the massive particles in the gauge theory, 

we will obtain just the logarithmic correction found in and reproduced in the 



above. If, however, we integrate out all the massive states above some characteristic 
mass scale (rricharact « jr), one can expect to have non-zero correlation functions 
between conformal operators of different conformal weight. This is also the structure 
of our mixed terms. Moreover, string loop corrections coming from the Born-Infeld 
action will be proportional to the above characteristic mass and will be negligible 
compared to the correction term for the correlators. 

5 Conclusions 

In the present paper, we have studied the Coulomb branch of SYM theory using 
AdS / CFT correspondence. Our considerations are based on the simplest example of 
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two stacks of N/2 parallel branes separated by a distance 2d which break SU{N) 
symmetry to S{U{N/2) x U{N/2)). Inspite of its apparent simplicity, the problem 
has turned out to be complicated but interesting. This complication mainly comes 
from the coupling between different n-wave dilatonic KK modes in the equations of 
motion. These equations are solved in the leading (1/A^) order. We have shown that 
the coupling between KK dilatonic modes causes the appearance of mixed correlation 
functions which are of order According to AdS/CFT correspondence principle, 
the operators correspond to tr{F'^X'^). It is obvous from the explicit form of 
the correlation functions that in the conformal limit {d 0, or Higgss vev tends to 
zero) the correction term in the diagonal correlators and the mixed correlation 
functions vanish. 

All the correlation functions have also deformation terms. The analysis in the 
phase space, after a Wick rotation, shows that these functions have a resonant 
structure - double poles at appropriate resonant frequencies. We interpret these 
double poles as the degenerate spectrum of glueballs. 

It will be interesting to study the correlation functions beyond this level of ap- 
proximation. We expect more conformal-like terms to appear in the correlations 
functions. We also have also to learn more about the spectrum of the glueballs and 
to look for a mechanism of deriving more realistic masses. 
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6 Appendix: The particular solutions of (|40])- 

In this Appendix we will find the particular solutions for the first order equations 
of motion. We have, in fact, to find particular solutions for three cases: 

^ 1 ^ / (n + 2?\ ^ I f^"^^^. 
fm + -MO-h+^—j^)fn{0 = { ein+2U) n = l,... (67) 
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Let us consider the following expression: 



V^H = H" + -M'{0- ( 1 + 



H{i) 



Since the right hand side of eq. ( p7| ) contain Bessel functions times powers of 
instructive to consider the function H of the form: 

After simple transformations and using the Bessel equation for we find: 



-a-2 



If we substitute for in the above equation: 
we find: 



a-2 



2ailm+i + ((a + mf - {n + 2f) /„ 



Let us choose the parameters a and m in the following way: 

a + m = n + 2 

Therefore: 

VnH = r-™+i2(n - m + 2)/^+i(0 
For m = n — 1 one find: 



d'' , 1 d f ^ {n + 2] 

de idi V e 



which solves the first case in eq. ([671). For m = n we have: 



d' ^ 1 d / ^ (n + 2)^ 
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To solve the second case one can use the properties of the Bessel function In+2- 

iln+2{i) = iUi) - 2{n + 2)/„+i(0 

and therefore: 

f 4+2(0 = f 4(0 - 2(n + 2)e7„+i(0 

Using the above we have the solution for the second case. 

To obtain the solution for the third case instead of eq. ( [7DD one can use the 
following property of Bessel functions: 

e4(0 = eX-i(0-"^/m(0 

Then we will have: 

Vr^H = [2aeX_i + ((a - mf -{n + 2f) /„] (75) 

Let us choose: 



With this choice we find: 
£ 1 d 



(72 + 2); 



e'/n+5(0 
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f/n+4(0 



(76) 



Now, consider the system of equations containing the modified Bessel functions 



{n + 2) 



(77) 



we find: 



rf2 1 d 

rf2 1 d 

df'^M 

d^ 1 d 

df^M 



- 1 



- 1 + 



(n + 2)2 



6 

4 



^X+i(0 

= fir„+4(0- 



(78) 
(79) 
(80) 
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